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ABSTRACT 

We present high signal-to-noise ratio Keck ESI spectra of the two quasars 
known to have Gunn-Peterson absorption troughs, SDSS J1030+0524 (z = 6.28) 
and SDSS J1148+5251 (z = 6.37). The Lya and Ly (3 troughs for SDSS 
J1030+0524 are very black and show no evidence for any emission over a red- 
shift interval of ~ 0.2 starting at z — 6. On the other hand, SDSS J1148+5251 
shows a number of emission peaks in the Ly (3 Gunn-Peterson trough along with 
a single weak peak in the Lya trough. The Lya emission has corresponding 
Ly (3 emission, suggesting that it is indeed a region of lower optical depth in the 
intergalactic medium at z = 6.08. 

The stronger Ly/3 peaks in the spectrum of SDSS J1148+5251 could conceiv- 
ably also be the result of "leaks" in the IGM, but we suggest that they are instead 
Lya emission from an intervening galaxy at z — 4.9. This hypothesis gains cre- 
dence from a strong complex of C IV absorption at the same redshift and from 
the detection of continuum emission in the Lya trough at the expected bright- 
ness. If this proposal is correct, the quasar light has probably been magnified 
through gravitational lensing by the intervening galaxy. The Stromgren sphere 
observed in the absorption spectrum of SDSS J1148+5251 is significantly smaller 
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than expected based on its brightness, which is consistent with the hypothesis 
that the quasar is lensed. 

If our argument for lensing is correct, the optical depths derived from the 
troughs of SDSS J1148+5251 are only lower limits (albeit still quite strong, with 
r(Lyo;) > 16 inferred from the Ly/3 trough.) The Ly/3 absorption trough of 
SDSS J 1030+0524 gives the single best measurement of the IGM transmission 
at z > 6, with an inferred optical depth r(Lyo;) > 22. 



1. Introduction 

The discovery of quasars at redshifts greater than 6 from the Sloan Digital Sky Sur- 
vey (SDSS; Fan et al. 2001, Fan et al. 2003) has opened a window for spectroscopic ex- 
ploration of the high redshift intergalactic medium. The detection of broad, black, Lya 
absorption troughs, as predicted by Gunn & Peterson (1965; hereafter GP), indicates that 
we are beginning to probe the era when hydrogen was reionized by an early generation 
of stars. Becker et al. (2001) detected the first complete GP trough in the spectrum of 
SDSSp J103027. 10+052455.0 (z = 6.28; hereafter SDSS J1030+0524). They compared the 
absorption to that seen in the spectra of lower redshift quasars and concluded that the frac- 
tion of neutral hydrogen in the IGM increases sharply at z ~ 6. Fan et al. (2003) report 
that the most distant known quasar, SDSS Jl 14816.64+525150.3 (z = 6.37 s ; hereafter SDSS 
J1148+5251), also shows a complete GP trough in its spectrum. Djorgovski, Castro, Stern 
& Mahabal (2001) argued that the spectra of even lower redshift quasars show indications 
of neutral regions in the IGM. 

There ensued a spate of theoretical papers discussing the details of how reionization 
takes place. But the most exciting new observational result is undoubtedly the discovery by 
the Wilkinson Microwave Anisotropy Probe (WMAP) that the microwave background has 
traversed a moderate optical depth (r = 0.17 ± 0.04) electron scattering medium, implying 
that the universe must have been ionized at z ~ 10-20 (Kogut et al. 2003). At first glance this 
appears inconsistent with reionization at z ~ 6, but since a relatively small neutral fraction 
(~ 1%) would suffice to produce black GP troughs, the results are not contradictory. There 
is mounting evidence that the ionization history of the intergalactic medium (IGM) was 
complex. Wyithe & Loeb (2003) suggested that the hydrogen in the IGM could have been 



8 This is a new redshift, different from that given in Fan et al. (2003), based on the spectrum presented 
in this paper. See §5.4 for a discussion of the z = 6.41 redshift derived by Willott, McLure & Jarvis (2003) 
for this object. 
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reionized twice, and Cen (2003) argued that two phases of reionization are likely under a 
wide range of conditions. In this scenario, the first reionization occurs at z ~ 15, driven by 
the formation of massive, zero-metallicity Population III stars (Venkatesan, Tumlinson & 
Shull 2003; Mackey, Bromm & Hernquist 2003). The ionizing radiation from massive stars is 
dramatically reduced, however, when Pop III supernovae seed the gas with metals (due to a 
reduction in the formation of massive stars, increased line-blanketing in hot star atmospheres, 
and increased mass loss rates from both hot stars and red giants.) Consequently the IGM 
recombines and remains relatively neutral, rim/nu > 0.1, until an increasing Pop II massive 
star population and a declining IGM density allow the second reionization to occur at z ~ 6. 
This neutral fraction is sufficient to produce the observed high optical depths in the Ly« 
GP troughs of high redshift quasars. 

To advance our understanding of the ionization history of the universe, we need high 
resolution, high signal-to- noise ratio spectra of the IGM absorption at z > 6. Recently 
several groups have discovered Ly a-emitting galaxies at z ~ 6.5 (Hu et al. 2002, Kodaira 
et al. 2003). These objects are interesting and relevant as direct probes of star and galaxy 
formation at z > 6, and the mere detectability of their Lya emission is an important clue 
to the ionization state of the IGM (Haiman 2002). However, their continua are far too faint 
(even redward of Lya) to make them useful for studies of the intervening IGM. Only high- 
redshift quasars offer the possibility of acquiring high quality information on the state of the 
IGM at high redshift. 

In this paper we present high quality spectra of the only two quasars known to have 
GP absorption troughs 9 , SDSS J1030+0524 (z = 6.28) and SDSS J1148+5251 (z = 6.37), 
taken with the ESI spectrograph on the Keck Telescope. The following sections describe 
the observations and data reduction (§2), discuss the analysis methods used (§3), present 
the properties of the IGM absorption in the spectra (§4), and discuss the results (§5). Ap- 
pendix A explores the effect of light travel time on the observed expansion of ionization 
fronts around high-redshift quasars. 

2. Observations 

Spectra of SDSS J1030+0524 and SDSS Jl 148+5251 were taken with the Keck Echel- 
lette Spectrograph and Imager (ESI; Sheinis et al. 2002) between January 2002 and February 

9 SDSS J1048+4637 at z = 6.23 (Fan et al. 2003) could also have a GP trough, but our recent spectra 
show that it has strong, intrinsic, broad absorption lines that make it less useful for studies of the IGM at 
z > 6 (Fan et al. 2003b). 
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2003. ESI covers the wavelength range 4000-11000 A with 10 Echelle orders having a con- 
stant resolution of R = A/AA ~ 10000 (~ 30 km/s). The data were taken with the slit at 
the parallactic angle in mostly good but variable seeing (from 0.6" to > 1", with mean seeing 
around 0.8"). The sky transparency ranged from photometric to light cirrus. A log of the 
observations is given in Table 1. 

The most challenging data reduction problem for these observations is subtracting the 
sky background, including the bright OH lines that dominate the near infrared atmospheric 
emission. The sky is considerably brighter than the object, and small errors in the sky 
subtraction can easily produce offsets in the emission that either create the appearance of 
light in the absorption trough or remove light that should be present there. To assist in the 
sky subtraction, we shifted the object along the slit for the various integrations, which helps 
to average over any small-scale flat-fielding, bias or dark current variations not corrected by 
the standard CCD calibrations. 

The sky brightness varies as a function of time in a complex way in the near infrared, 
with the atmospheric OH emission varying rapidly and independently of other sky emission. 
In addition the ESI spectrograph, though a highly stable instrument, does undergo small 
changes over time, with the wavelengths shifting slightly from one observation to the next. 
Presumably these changes (which are typically at the level of a few tenths of a pixel or less) 
are due to minor thermal and gravity-induced flexures. 

We extract a sky-subtracted spectrum from each independent observation. The separate 
spectra are then combined using weights to optimize the signal-to- noise ratio of the result. It 
is necessary to be very careful with the weights for combining to avoid any bias in the results. 
A simple unweighted average is safely free of bias but leads to a much noisier spectrum than 
an optimally weighted sum. On the other hand, using weights l/o" 2 with a derived from 
the noise estimates of the individual spectra leads to a very noticeable bias toward lower 
flux values, because the estimated noise is larger for brighter pixels. In the black absorption 
troughs this leads to fluxes that are systematically negative. 

The dangers of using a weighted combination for Poisson data are well-documented (for 
a nice analysis, see Wheaton et al. 1995). The cure for the bias is to choose weights that are 
completely independent of the data. We have combined the spectra using weights that are 
derived from the sky noise in adjacent pixels (i.e., along the slit). Note that using the sky 
noise from the same pixel would still lead to a small bias, because pixels that have positive sky 
fluctuations (making the residual counts a little low) would be given slightly lower weights, 
thus leading to a positive bias in the sum. Using the sky noise in neighboring pixels produces 
a completely unbiased result that has a nearly optimal signal-to-noise ratio. An overall scale 
factor is determined for each spectrum to match the fluxes, thereby correcting for variable 



- 5 - 



atmospheric transparency. 

Our extraction algorithm is highly customized for ESI data, which have a number of un- 
usual characteristics (e.g., highly curved orders, a strongly non-linear wavelength dispersion, 
and some odd notches in the sensitivity in the bluer orders.) The images are first geomet- 
rically rectified to produce orders with orthogonal wavelength and spatial coordinates lying 
along rows and columns. This removes the very large curvature of the ESI Echelle orders 
and also corrects small variations in the wavelength dispersion along the slit that manifest 
themselves as tilted or curved sky lines. The order curvature removal also automatically 
corrects the spectral trace for differential atmospheric refraction using the airmass of the 
observation; this approach is more effective than tracing the orders directly because high- 
redshift quasar spectra have long stretches with no detectable light, foiling standard tracing 
algorithms. 

The resampling algorithm used in the rectification is a flux-conserving interpolation 
scheme that is accurate to second order in the pixel shifts (as opposed to the more com- 
monly used schemes that are only first order, making them considerably more dispersive.) 
The interpolation method is based on the high-order, monotonic, flux-conserving advection 
schemes that are widely used in modern hydrodynamic simulations (van Leer 1977). 

Extraction of spectra from the rectified orders is straightforward. We use an optimal 
extraction algorithm with a second-order fit to the sky along the slit and a spatial profile that 
is assumed constant for each order. For orders where the profile is not detected, the trace 
position and width are predicted based on those orders that are detected. The wavelength 
scale derived from calibration lamps is adjusted using the positions of sky lines for each 
observation. The final spectra for the quasars are shown in Figure 1. 

3. Analysis of the Spectra 

High signal-to-noise ratio spectra promise two improvements in our understanding of the 
GP troughs. By reducing the noise in the black regions of the spectrum, they improve our 
lower limits on the absorbing optical depth and may allow us to detect faint light in the dark 
troughs. But since the optical depth limits increase only as the logarithm of the signal-to- 
noise ratio, which itself increases only as the square root of the integration time, it is difficult 
to make dramatic improvements in the limits. The main benefit of high signal-to-noise ratio 
spectra is that they allow us to set strong detection limits over narrower bands in the spectra. 
The optical depth distribution is very unlikely to be a smooth function of redshift; instead, 
at the transition between a neutral and an ionized IGM, there will be ionized, partially 
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Fig. 1. — Keck ESI spectra of the two highest redshift quasars. The spectra have been binned 
to a resolution R = AA/A = 2500 (120 km/s). The Gunn- Peterson absorption troughs of 
Lya and Ly (3 are indicated, extending from z = 5.975 to 6.165 (SDSS J 1030+0524) and 
z = 6.075 to 6.320 (SDSS J1148+5251). 
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transparent bubbles sprinkled along the line of sight. To detect these bubbles we need to be 
able to recognize narrow emission spikes in the trough where the quasar light leaks through 
the IGM. The new spectra presented in the paper are far better for that purpose than our 
earlier observations (which were 30 min exposures taken with ESI.) 

We begin with an approach similar to that used by Becker et al. (2001). The residual 
fluxes in the GP troughs of Lya and Ly/3 are measured directly from the spectra; the 
fluxes are weighted using the sky noise from adjacent pixels (as described above for the 
extraction) to improve the signal-to-noise ratio while avoiding bias in the sums. Then the 
mean transmission and optical depth of the IGM are determined by comparing the residual 
flux to the (unabsorbed) quasar continuum at the positions of the troughs. 

Because very little neutral hydrogen is required to create a GP trough in Lya (a 1% 
neutral fraction will suffice), the presence of a GP trough due to Lya provides only a lower 
limit to the amount of neutral hydrogen. The GP trough from the Ly/3 transition can be 
a more sensitive test because the Ly (3 cross-section is a factor of 5.27 lower than Lya. 
Interpretation of the Ly/5 trough is complicated by the fact that it is overlaid by the Lya 
forest at z ~ 5, which is unpredictable along a given line of sight due to cosmic variance. In 
spite of that, using fairly secure assumptions, the Ly/3 trough can give much stronger limits 
on the neutral hydrogen fraction. We consequently measure residual fluxes in both the Ly a 
and Ly (3 absorption regions. 

We also considered using the Ly7 absorption since Ly7 has an even lower intrinsic 
optical depth. However, the Ly7 GP trough suffers from overlying absorption by both Lya 
and Ly (5 at lower redshifts and consequently does not give additional useful information. 

It is necessary to estimate the brightness of the quasar continuum at the Lya and 
Ly/? troughs in order to measure the transmission of the IGM. For this purpose we use 
the far UV quasar power law spectrum of Telfer, Zheng, Kriss, & Davidsen (2002) and the 
composite spectrum from the Large Bright Quasar Survey (LBQS; Francis et al. 1991 10 ). 
The templates are matched to the observed spectra at a rest wavelength of ~ 1290 A. The 
Telfer et al. power law has a break near 1280 A and rises less steeply into the UV than a 
simple extrapolation of the continuum redward of Lya. In the near UV (A > 1280 A), the 
spectral slope is a^uv = —0.72 (f u oc v Q ), while in the far UV (A < 1280 A), 

&EUV — — 57 

for radio-quiet objects and —1.96 for radio-loud objects. While there is no evidence that 
either of our quasars is radio-loud, we compute continuum levels using both EUV slopes in 



10 We use the improved LBQS template created by S. Morris, described in Brotherton et al. (2001) and 
available at http://sundog.stsci.edu/first/QSOComposites. Note that the LBQS composite is very similar 
to the FIRST Bright Quasar Survey composite derived by Brotherton ct al. 
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order to represent the range of continua seen in quasars. 

Broad emission lines can also contribute to the continuum. Our spectra are displayed 
in Figures 2 and 3 along with the LBQS composite spectrum. For these figures, our spectra 
have been denoised using a wavelet shrinkage algorithm (Donoho & Johnstone 1994) that is 
optimized for the case where the noise amplitude is variable but known (heteroscedastic data, 
in the parlance of statisticians). We use a Haar (1910) wavelet transform, which allows error 
propagation and has the great benefit of keeping the noise in transform space uncorrelated at 
each scale before the filtering is applied. Cycle-spinning (Coifman & Donoho 1995) is utilized 
to make the result shift-invariant and to remove the artifacts generated by the square-wave 
Haar basis functions. 

Denoising thus uses a variable smoothing length to suppress the noise as much as possible 
on all scales while remaining consistent with the spectrum given the estimated noise array. 
The resulting spectrum is heavily smoothed in flat regions (such as the dark troughs), which 
suppresses the noise. (Values in the smoothed areas are consequently highly correlated.) 
The smoothing threshold is set to 3a, so only structures more significant than that remain 
in the spectrum. Note that high resolution structures are preserved (compare the absorption 
on the blue edge of the Lja emission with Fig. 1) while the noise in both the continuum 
and the troughs is dramatically reduced. 

For SDSS J1030+0524 (Fig. 2), the LBQS composite is clearly a good match to the 
spectrum, so we estimate the continuum in the Ly a and Ly (5 absorption troughs using the 
mean flux in the relevant wavelength window from the LBQS spectrum. On the other hand, 
SDSS Jl 148+5251 shows considerably weaker emission lines than the composite (Fig. 3); we 
therefore use the Telfer et al. power law continuum alone as the estimate of the brightness 
at the absorption troughs. The continua estimated by using these different templates differ 
by about 30%. 

4. Results 

4.1. SDSS J1030+0524 

For SDSS J1030+0524, our results can be described very simply: the GP absorption 
troughs of both Ly a and Ly (5 are very black. There is no evidence of any residual light in 
the troughs to the limit of our observations. The Lya and Ly/3 absorption trough details 
are shown in Figure 4. 

The results are summarized in Table 2. The band limits were chosen to exclude the first 
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Table 1. Keck ESI Observation Log 



OUUI Ct: 




Dill VV 1QL11 


in uiiiuei ui 






yyyy/ mm/ dd 


(") 


Exposures 


(hr) 


SDSS J1030+0524 


2002/01/10 


1.00 


6 


2.00 


SDSS J1030+0524 


2002/01/11 


1.00 


5 


1.67 


SDSS J1030+0524 


2002/01/12 


1.00 


12 


4.00 


SDSS J1030+0524 


2002/01/13 


0.75 


2 


1.50 


SDSS J1030+0524 


2002/02/07 


0.75 


2 


1.50 


SDSS J1148+5251 


2002/12/04 


1.00 


9 


3.00 


SDSS J1148+5251 


2002/12/05 


1.00 


10 


3.33 


SDSS J1148+5251 


2003/02/04 


1.00 


14 


4.67 



Table 2. Absorption Properties of z > 6 Quasars 



Object 


GP Trough 








Zmax 


Residual Flux 


Continuum 


IGM Transmission 


Lya Optical 






(A) 


(A) 






(10~ 20 ergs~ 


'cm- 2 !- 1 ) 


Fraction 


Depth 


SDSS J1030+0524 


Ly a 


8510 


8710 


6.00 


6.17 


1.6 ± 1.4 


1570 


0.0010 ± 0.0009 


> 6.3 a 


Zem = 6.28 


Ly/3 


7180 


7349 


6.00 


6.17 


0.8 ± 1.7 


200 b 


0.0043 ± 0.0088 


> 22. 8 a 


SDSS J1148+5251 


Lya 


8630 


8900 


6.10 


6.32 


3.6 ± 0.8 


1730 


0.0021 ± 0.0005 


6.2 ±0.2 


^em — 6.3T 


Ly/3 


7282 


7509 


6.10 


6.32 


9.6 ± 1.0 


210 c 


0.0465 ± 0.0047 


16.2 ±0.6 




Ly/3 


7282 


7509 


6.10 


6.32 


9.6 ± 1.0 


190 d 


0.0510 ± 0.0052 


15.7 ±0.6 




Lya 


8510 


8630 


6.00 


6.10 


10.9 ± 1.1 


1740 


0.0063 ± 0.0007 


5.1±0.1 




Ly/3 


7180 


7282 


6.00 


6.10 


63.4 ±1.4 


210 c 


0.3020 ± 0.0067 


6.3 ±0.1 




Ly/3 


7180 


7282 


6.00 


6.10 


63.4 ±1.4 


190 d 


0.3350 ± 0.0074 


5.8 ±0.1 



a lcr lower limit to the optical depth. 

b Derived from the LBQS composite spectrum with a Lya forest absorption factor of 0.12. 

c Derived from the Telfer et al. (2002) radio-quiet spectral index, oleuv = —1-57 with a Lya forest absorption factor of 0.11. 
d Derived from the Telfer et al. (2002) radio-loud spectral index, oleuv = —1-96 with a Lya forest absorption factor of 0.11. 
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Fig. 2. — Denoised, full-resolution spectrum of SDSS J1030+0524 with matched templates 
from the LBQS and Telfer et al. (2002). The denoising algorithm is described in the text. 
The template is a very good match to the quasar redward of the Ly a IGM absorption. 
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Fig. 3. — Denoised, full-resolution spectrum of SDSS Jl 148+5251 with matched templates 
from the LBQS and Telfer et al. The hydrogen emission lines are much weaker than those 
in the template. Our best redshift determination for this object, z = 6.37, comes from 
matching the features redward of 9500 A with the template. 
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Fig. 4. — Closeup of the Lya and Ly/3 GP troughs, aligned in redshift, at a resolution of 
R = 2500 (120 km/s). Note that these spectra have not been denoised, so the pixel values 
are independent. The top axis shows the wavelength for Lya, and the shaded band shows la 
errors. The absorption is extraordinarily black for SDSS J1030+0524, but SDSS J1148+5251 
shows evidence for weak emission in the trough, especially at Ly (5. The strong peaks in the 
Ly (3 spectrum at z — 6.02 and 6.06 are discussed in the text; they could be Lya emission 
from a foreground galaxy at z — 5. 
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bright peaks on either end of the Ly (3 GP trough and are identical in redshift for the two 
troughs. The residual fluxes for both the Ly a and Ly j3 troughs are consistent with zero. 

For the Ly/3 trough, the continuum estimate includes a reduction factor of 0.12 for 
the overlying Lya forest absorption at z = 5 (as determined by Becker et al. 2001) The 
r(Ly P) values are converted to equivalent Lya optical depths by multiplying by 5.27, the 
ratio of the oscillator strengths of the two lines. Thus the Ly (3 GP trough measurements for 
SDSS J1030+0524 imply r(Lya) > 22.8 (lcr). The errors quoted in the table reflect only 
the statistical uncertainties. As found by Becker et al., the strongest limit on the optical 
depth comes from the Ly f3 trough, despite the additional uncertainty that comes from the 
overlying Lya absorption. 

These results are consistent with the measurements of Becker et al. (2001) (confirmed 
by Pentericci et al. 2002) but are naturally stronger given the much improved signal-to-noise 
ratio of the data. Our limits on the residual flux in the GP troughs and the Lya optical 
depth of the IGM are lower by the expected factors given the longer integrations reported 
here. The limits remain consistent with zero flux in the GP troughs of SDSS J1030+0524. 

There is no evidence for any significant peaks in either trough. Figure 5 displays the 
signal-to-noise ratio of the spectra in the troughs. The distribution is consistent with zero- 
mean noise in the troughs; the reduced xl f°r the combined Lya and Ly/3 troughs is 0.97. 

4.2. SDSS J1148+5251 

The GP troughs for SDSS Jl 148+5251 are also shown in Figures 4 and 5. The mean 
transmission properties of the troughs are given in Table 2, which includes the transmission 
measured for both a high redshift range (6.1 < z < 6.32) and a lower redshift interval 
(6.0 < z < 6.1). Note that we now quote a measurement, not a lower limit, on r, as flux has 
been detected in the troughs. The correction factor for Ly a forest absorption overlying the 
Ly (3 GP trough was taken to be 0.11, from Figure 2 of Becker et al. (2001). The la errors 
on the transmission are smaller than for SDSS J1030+0524, both because the seeing and 
photometric conditions were slightly better for SDSS J1148+5251 and because its continuum 
is brighter. 

In contrast to SDSS J1030+0524, there is light detected in the GP troughs of SDSS 
Jl 148+5251 both in narrow spikes and in broader bands. The mean transmission nonetheless 
remains low. Figure 6 summarizes the Lya optical depth as a function of redshift, including 
additional high-redshift quasars from Fan et al. (2001). There is a noticeable steepening of 
the t(z) relationship at z ~ 5.8. 
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Fig. 5. — Signal-to-noise ratios for the Ly a and Ly (3 absorption troughs in SDSS J1030+0524 
and SDSS J1148+5251. 
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Fig. 6. — Lya optical depth as a function of redshift from high-redshift quasars. Error bars 
(and the lower limit for SDSS J1030+0524) are 2a. Data for SDSS J1030+0524 and SDSS 
J1148+5251 are from this paper; other data are from Fan et al. (2001). The line shows the 
expected evolution when the density of Ly a forest clouds scales as (1 + zfl 2 (also from Fan 
et al.) The implied r(Lya) values from the Ly/3 GP trough transmission are substantially 
stronger: r > 22 (SDSS J1030+0524) and r = 16 (SDSS J1148+5251). We believe that 
highest redshift bins for SDSS Jl 148+5251 are probably contaminated by light from an 
intervening galaxy and so should be considered lower limits. 
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A closer examination of the troughs reveals some interesting features. The Ly a trough 
extends from z = 5.93 to 6.32 and looks dark at high resolution except for a weak emission 
spike at z — 6.08. The reality of that peak is supported by its coincidence with a Ly/3 
emission peak at the same redshift (Fig. 4). We believe this is likely to be a real detection 
of an ionized bubble in the IGM that permits quasar light to leak through at the redshift of 
the bubble; this feature is discussed in detail below (§5.1). 

Several emission spikes are detected in the Ly/3 trough as well. Two especially strong 
peaks are seen at 7205 A and 7236 A (Ly/3 redshifts of z = 6.03 and 6.06). There also 
is significant non-zero flux in the troughs even between the peaks (best seen in Fig. 5). 
If this flux is quasar light transmitted by the IGM, it certainly suggests strongly that the 
IGM is not very neutral at z ~ 6. That interpretation is not so clear cut, however. We 
argue below (§5.2) for the alternative interpretation that both the Lya and Ly/3 troughs 
are contaminated by light from an intervening galaxy at z — 5. If that suggestion is correct, 
the optical depths derived from the troughs of SDSS J1148+5251 are only lower limits. 



5. Discussion 

5.1. The z = 6.08 emission feature in SDSS J1148+5251 

The universe at z ~ 6 was very likely ionized by galaxies rather than by quasars, since 
there do not appear to be enough quasars to do the job (Fan et al. 2001). There must have 
been regions in the preionized universe where a concentration of unusually luminous galaxies 
(possibly the precursor to a massive cluster) was capable of ionizing the nearby IGM, creating 
an ionized bubble. It is therefore not surprising to find weak emission spikes embedded in 
the mostly black GP troughs (Figs. 4 and 5) where the line-of-sight to the quasar passes 
through small "holes" in the IGM Ly a and Ly f3 transmission created by ionized bubbles. 
That can also explain the detection of Lya-emitting galaxies at z ~ 6.5 (Hu et al. 2002, 
Kodaira et al. 2003), which in this picture ought to be found in concentrations with other 
galaxies. 

According to the scenario presented by Cen (2003), the universe is first ionized at z ~ 15 
by Pop III stars and then recombines to become substantially neutral, with neutral fractions 
remaining above 0.1 until the second reionization occurs at z ~ 6. At z = 6 to 7 the neutral 
fraction 71hi/?ih is typically ~ 0.15. As Cen points out, that makes it easier to detect high 
redshift galaxies; it also makes it easier for lower luminosity sources to create ionized bubbles. 

The peak seen in the Lya trough of SDSS Jl 148+5251 at 8600 A (z = 6.08) has an 
intensity about 8% of the extrapolated (pre-absorption) continuum, implying an IGM optical 
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depth of r(Lya) = 2.5. The optical depth of this window at Ly/3 would be only about 0.5, 
which would produce a very strong peak in the Ly/5 trough; but since the overlying Lya 
forest absorption is expected to reduce that by a factor of 0.11 (corresponding to an overlying 
optical depth r = 2.2), the expected amplitude of the Ly (3 trough peak is in fact very similar 
to that observed, albeit with a rather large uncertainty in the Ly a forest attenuation. 

How large would an ionized region need to be to produce this feature? The damping 
wings from neutral hydrogen on either side of the ionized bubble create a high optical depth 
through the region unless the bubble is relatively large. This led Barkana (2002) to question 
Djorgovski et al.'s (2001) argument that dark regions in the spectrum of a lower redshift 
were created by neutral regions of the IGM. Miralda-Escude (1998; see also Barkana 2002) 
derived the Lja damping wing optical depth for hydrogen extending from z\ to z 2 when 
observed at a wavelength A a (l + z): 
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Given an ionized region at redshift z, Eqn. (1) allows us compute the combined damping 
wing optical depth from the neutral regions on either side of the bubble. Although the 
constant factor multiplying r G p in Eqn. (1) is small, the leading term in I(x) can be very 
large near z = Z\ or z 2 , so r damp can be substantial. 

Figure 7 displays the resulting optical depth for a bubble located at z — 6.08. The 
neutral IGM is assumed to extend from z = 6 to z = 6.32, but the results are relatively 
insensitive to those limits. The bubble diameter required to produce a window with r = 2.5 
is 1.33 Mpc if the IGM is completely neutral outside the bubble. If the IGM is already partly 
ionized, the required bubble size is much smaller (0.22 Mpc if fiHi/% = 0.15). 

Barkana (2002) computed the diameter of an H II bubble as a function of the halo 
mass M of the ionizing source (including the escape probability for ionizing photons in a low 
metallicity galaxy) to be 

Bs = 0,50Mpof^f^J^V /3 , (4) 
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Fig. 7. — Optical depth from damping wings of Lya as a function of diameter for an 
ionized bubble at z — 6.08. For small bubbles in a completely neutral IGM (solid line), 
the damping wings strongly suppress the transmission. If the IGM is already partly ionized 
(n H i/n H = 0.15; dashed line) then the optical depths are smaller by the same factor. The 
dots show the bubble diameter required to produce the observed peak in the Ly a trough of 
SDSS J1148+5251 (1.33 Mpc for a neutral IGM or 0.22 Mpc for a partially ionized IGM.) 
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where A^ on is the number of ionizing photons per baryon that escape the galaxy. We can 
invert this to get the mass required to produce a bubble of a given diameter: 

We conclude that a plausible star-forming galaxy (or group of galaxies) along the line of 
sight can produce the observed peak at z = 6.08 in the Lya and Ly/5 GP troughs of SDSS 
J1148+5251, even if the surrounding IGM is completely neutral. Obviously if the IGM is 
partly ionized, the required bubble diameter would be much smaller (Fig. 7) and so the 
required galaxy mass would also be far smaller (< 1O 9 M ). With future instrumentation it 
will be interesting to search the vicinity of the quasar for evidence of star-forming galaxies 
at z — 6.08 that are the source of the ionizing continuum radiation required to make this 
bubble. 



5.2. Evidence for an intervening z = 5 galaxy 

The strong spikes seen in the Ly/3 trough of SDSS Jl 148+5251 at 7205 A and 7236 A 
are more difficult to explain as the result of ionizing bubbles in the IGM. The strength of 
these peaks is not easily reconciled with the absence of corresponding Lja peaks, given the 
expected overlying Lya forest absorption. The peak at 7205 A has an amplitude 8.7 ± 0.3 
(xlO -18 ergs _1 cm~ 2 A ), while the brightest plausible continuum at that wavelength is 24 
(from the LBQS template spectrum in Fig. 3). If we assume that the 7205 A peak is quasar 
light leaking through the z=6 Ly (5 and z=5 Ly a absorption, the implied combined optical 
depth is t(Lj a@z=5) + r(Ly (3@z—6) = 1.0. But there is no hint of a Lya peak at 8539 A: 
the flux limit there is 0.1 ± 0.2 with an estimated continuum level of 17, implying 3a upper 
limits r(Lya@z=6) > 3.3 and r(Ly /3@z=6) = r(Lya@z=6)/5.27 > 0.6. The overlying 
Lya forest absorption must then be t(Lj a@z=5) < 0.4. 

We conclude that if the 7205 A peak is the result of an ionized bubble at z — 6.03, it 
must fall in a highly transparent window in the z = 5 Ly a forest that barely attenuates the 
continuum at all, allowing > 65% of the light through. High resolution spectra of the Lja 
forest at z ~ 5 (Songaila, Hu, Cowie & McMahon 1999; Djorgovski et al. 2001; Songaila & 
Cowie 2002) do show small transparent windows, but they are rare. The spectra presented 
in this paper are typical. There is one transparent window seen at 7800 A(z = 5.4) m our 
spectrum of SDSS J1030+0524 (Fig. 2), and no windows at all are seen in SDSS Jl 148+5251. 
In the 5 < z < 5.5 Lya forest, only 0.5% of the bandpass has optical depths less than 0.4. 
Songaila & Cowie (2002) discuss the distribution of the transmitted fraction as a function 
of redshift; using their analytical model, at z = 5 only 1.7% of the spectrum is expected to 
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have a transmission > 0.65. This fraction drops to only 0.6% if we use the la upper limit, 
r(Lya@z=5) < 0.2. 

We conclude that it is problematic to reconcile the strong peak at z ~ 6.03 in the Ly/3 
trough with the absence of a corresponding Lya peak. Unless the Lya forest is extraor- 
dinarily transparent at the position of the peak, we would expect a strong, easily detected 
peak in the Lya GP trough. 

There is, however, another possible explanation for the strong peaks in the Ly (5 trough 
of SDSS J1148+5251. A careful examination of the quasar spectrum (Fig. 8) reveals a strong 
C IV absorption feature at 9200A (z = 4.943). That suggests a second interpretation for the 
spikes, namely that they are Lya emission from a foreground galaxy at z = 4.943. Usually 
one does not see emission associated with C IV absorbers; but ordinarily such emission lines 
would be lost in the Lya forest and not noticed. This particular emission line, falling as it 
does in the Ly/5 trough, stands out like a sore thumb. Moreover there are multiple C IV 
absorption systems and multiple (putative) Lya emission components seen, suggesting the 
presence of a fairly massive mass concentration. 

Is it plausible that the emission in the Ly/3 trough could be from an intervening z = 4.9 
galaxy? The implied Ly a luminosities are similar to those of Ly a emitters detected in the 
Subaru deep field at z — 4.86 (Ouchi et al. 2003). Moreover, the velocity offset between the 
Lya emission and the C IV absorption seen in Fig. 8 is in the range typically observed for 
systems that show both Lya emission and metal line absorption (~ 500 km/s, Shapley et 
al. 2002). 

The main argument against the intervening emission model is that Lya emitters are 
sufficiently rare on the sky that it is unlikely to find one along a random line-of-sight. But 
perhaps this is not a random line-of-sight. Figure 8 shows that there are both multiple Lya 
emission components and multiple C IV absorption components present in the spectrum 
(note particularly the structure in both components of the strong C IV doublet at 9200 A.) 
The two presumed Lya lines are split by 1200 km/s, and the C IV structure extends over at 
least 250 km/s. It is possible that there is a massive structure (a protocluster?) along the 
line-of-sight that is amplifying the quasar through gravitational lensing. If that is correct, 
then the likelihood of finding a Lya emitter in front of the z = 6.37 quasar is in fact not 
so small, since the foreground object assisted in the discovery of the quasar by boosting its 
apparent brightness. Even lensing by a factor of two might create a strong discovery bias 
since the luminosity function is steep for bright quasars (Wyithe & Loeb 2002; Comerford, 
Haiman & Schaye 2002). 

For the more familiar case where a nearby (z < 1) galaxy is lensing a high-redshift 
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Fig. 8. — Comparison of emission in the Ly/3 GP trough of SDSS J1148+5251 with absorp- 
tion by the C IV 1548.20, 1550.77 doublet, under the assumption that the emission lines are 
Lja from an intervening galaxy at the same redshift as the C IV absorption, z = 4.943. 
Probable C IV absorption systems (having both the correct separation and the correct opti- 
cal depth ratio) are marked by brackets in the lower panel. The panels are aligned in redshift 
so that the Ly a peak would fall midway between the C IV doublet components for perfect 
alignment. The velocity axis is registered to the strongest C IV at z = 4.943. The data are 
shown at the full resolution of ESI (a substantially higher resolution than is shown in Fig. 4.) 
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quasar, the lensing magnification is limited to less than ~ 2 because stronger lensing produces 
multiple images with a noticeably wide splitting. But when the lens redshift is z ~ 5, as 
in the scenario proposed here, the Einstein ring radius is much smaller and so the splitting 
could not have been detected from the ground (in agreement with the imaging observations 
by Fan et al. 2003.) Figure 9 shows the magnification and observed image positions as a 
function of impact parameter for a singular isothermal sphere galaxy model with a velocity 
dispersion of 250 km/s. This velocity dispersion is consistent with the velocity structure 
seen in the strong C IV absorption system and is slightly larger than the dispersion expected 
for an galaxy (v = 220 km/s). In this geometry, the multiple image separation is only 
0.3 arcsec, which would be difficult to detect even if the components were of equal brightness. 

Note also that since there are multiple C IV absorption and Ly a emission components, 
it is possible that there are multiple lensing galaxies along the line-of-sight to this source. 
That could possibly increase the amplification factor further, although it would also lead to 
larger splitting and shifts. We coadded the two-dimensional spectra of SDSS J1148+5251 
for the Echelle order that includes the Ly f3 trough to search for evidence that the peaks are 
either broader than the quasar trace or are shifted relative to the center of the trace. Neither 
effect was seen. The width of the peaks agrees with the width of the rest of the order, and 
the central positions agree to better than 0.1" (which is the approximate uncertainty in the 
trace position.) 

This appears to be the most problematic aspect of interpreting these peaks as Lya 
emission from a foreground lensing galaxy: we see no direct evidence any any images (in- 
cluding the spectrum) for lensing. That does not contradict the lensing hypothesis given the 
small lens scale at z — 5, but neither does it provide any additional support. And there is 
certainly reason to think it unlikely that we would stumble across a lens with such a small 
Einstein ring radius, since the area on the sky that is lensed by such objects is rather small. 
But as we discuss below, the additional supporting evidence in favor of the hypothesis is 
sufficient to make a reasonably strong case. 

5.3. Continuum emission in the troughs of SDSS J1148+5251 

A puzzling aspect of the Lya and Ly (3 GP trough transmissions for SDSS J1148+5251 
is that the Ly a residual flux is far higher than expected from the Ly (3 flux. This is reflected 
in Table 2 by the fact that r(Lya) inferred from the Ly/3 trough is far higher than the 
optical depth actually measured in Lya. Either there is too much light in Lya trough or 
there is too little in the Ly (3 trough (which presumably would require a great deal of extra 
absorption in the z ~ 5 Lya forest.) Neither explanation appears particularly palatable 
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Fig. 9. — Gravitational lensing magnification (top) and splitting (bottom) of multiple com- 
ponents for a singular isothermal sphere with a velocity dispersion of 250 km/s. The solid 
line shows the primary, brighter component; the fainter second image that appears for impact 
parameters less than the Einstein radius is shown with a dashed line. The total magnification 
when the two components are combined is shown with a dotted line. Note that the Einstein 
ring radius is quite small for a quasar at z — 6.4 being lensed by a galaxy at z — 5.0, so 
even high magnification factors are not ruled out by the absence of multiple components in 
the image. 
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based on a conventional model of a quasar with IGM absorption. 

An intervening Ly a-emitting galaxy at z = 5 offers a more natural explanation: the Ly a 
and Ly (3 troughs are contaminated by continuum emission from the galaxy. We find that the 
observed brightness of the continuum is fully consistent with the hypothesis of an intervening 
galaxy, which strongly supports our suggestion. Wyithe & Loeb (2002) pointed out that a 
lensing galaxy could contaminate the black GP trough, although they were concerned mainly 
with lensing by low redshift galaxies (which is a priori considerably more likely.) 

Between 8420 A and 8900 A (z — 5.93 to 6.32), the only obvious emission in the Ly« 
trough is the peak at z = 6.08 discussed above (Fig. 4). If we sum the light in that wavelength 
range, excluding the peak (8600 < A < 8620), we get an estimate for the continuum level 
of 3.9 ± 0.6 x 10~ 20 ergs _1 cm~ 2 A . We take that as an estimate of the brightness of the 
galaxy continuum. The AB magnitude of the flux detected in the Lya trough is 26.4, which 
corresponds to an absolute AB magnitude of —21.9 at a rest wavelength of 1450 A if the 
light comes from a galaxy at z — 5. 

The observed equivalent width (EW) of the peak at 7205 A is then 930 A, and the 
rest-frame EW is 160. This EW is perfectly consistent with the EW distribution determined 
by Malhotra & Rhoads (2002) for Ly a-emitting galaxies at z — 4.5. It is slightly lower than 
the median EW for their sample, as is expected for a more massive galaxy (Malhotra & 
Rhoads, private communication). 

5.4. The extent of the Stromgren sphere 

Many authors have discussed the possibility that high-redshift quasars are selected be- 
cause they are strongly lensed. Wyithe & Loeb (2002) and Comerford, Haiman, & Schaye 
(2002) determined how lensing might bias the quasar luminosity function. Haiman & Cen 
(2002; hereafter HC) considered the specific case of lensing of SDSS J1030+0524; they argue 
that if the quasar is embedded in a largely neutral IGM, the large extent of its H II region 
implies that it must in fact have a large ionizing luminosity and so cannot have been too 
strongly lensed. If the ionizing luminosity were too low, the IGM near the quasar would 
have remained neutral and so the Lja emission line would have been much more heavily 
absorbed on the blue side. This argument also sets a minimum age for the quasar, since it 
must have been in existence long enough to ionize a substantial Stromgren sphere. 

The HC result is directly relevant to the question of whether SDSS Jl 148+5251 could be 
strongly gravitationally lensed. On the blue side of Ly a, emission is detectable to z — 6.32 
(Fig. 4), a comoving distance of 2.8 Mpc from the quasar. (Note that this distance is 
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somewhat uncertain due to the difficulty in determining an accurate redshift for the quasar; 
this is discussed further below.) According to HC — see also Madau & Rees (2000), Cen 
& Haiman (2000) and references therein — the radius of a Stromgren sphere ionized by a 
high-redshift quasar as a function of time is approximately 

(\ x / 3 -l 
— tt— j — — ] , (6) 

lO^s" 1 10 7 yr J V 7 - 37 / 

where N p h is the ionizing photon luminosity, £q is the age of the quasar, and as above we use 
cosmological parameters Hq = 70, Q\ = 0.7, Q m = 0.3, and = 0.04. This formula results 
from matching the total number of ionizing photons emitted to the number of hydrogen 
atoms within a sphere, and it ignores both recombinations (which reduce the radius) and 
the Hubble expansion (which increases it.) It uses the mean IGM density assuming that the 
IGM contains nearly all the baryons; if the universe is locally overdense by a factor D, the 
radius is reduced by D^ 1 ^. 

Inverting Eqn. (6) to solve for the quasar parameters, we find 

X "" tQ = 0.063 ( Rs 1 + ZQ X . (7) 
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The ionizing luminosity for SDSS J1148+5251 (assuming no lensing) is rather uncertain. 
We have estimated it using two different quasar spectral templates. The Telfer et al. (2002) 
quasar template, which we have used for the Ly a and Ly (5 continuum normalization, gives 
an ionizing photon luminosity of N p h = 0.9-1.3 x 10 58 s _1 , where the smaller value assumes 
the EUV spectral index for radio-loud quasars (oiEuv — —1-96) and the larger assumes the 
radio-quiet quasar index (oceuv — —1-57). The Elvis et al. (1994) template, which was used 
by HC, gives the substantially smaller value N p h = 0.2 x 10 58 . We think the Telfer et al. 
spectrum is preferable but report results using both spectra below. The corresponding range 
of values for t Q is 5-7 x 10 5 yrs or 3x 10 6 yrs (for the Telfer and Elvis templates, respectively.) 

Such a short lifetime for SDSS J1148+5251 would be surprising, since its high luminosity 
implies a large black-hole mass which should take a considerable length of time to assemble. 
The e-folding timescale for an accreting black hole is 

t acc = Ax 10 7 yr(^ r )^ 1 , (8) 

where e is the radiative efficiency and rj is the ratio of the quasar's luminosity to the Ed- 
dington luminosity (Haiman & Loeb 2001). This long natural timescale makes it unlikely 
that we will observe quasars with ages as short as a few hundred thousand years unless the 
radiative efficiency e is very low, allowing the quasar to accrete quietly and rapidly. 
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On the other hand, following HC's argument, if the quasar is lensed then its true lumi- 
nosity could be substantially smaller. With a plausible age of tq ~ 10 7 yrs, SDSS J1148+5251 
could be easily lensed by a factor of 3 or more and still be able to ionize the IGM out to the 
observed distance of 2.8 Mpc. 

We note that for the sizes and ages derived for this object, the Stromgren sphere is still 
in its early phase of rapid relativistic expansion. Even in that stage of evolution, however, 
Eqn. (6) accurately describes the observed size of the H II region as measured by absorption 
along the line-of-sight. See the Appendix for further details. 

One more point that is worth making is that these discrepancies become even greater if 
the IGM is not completely neutral at z ~ 6. The expected value of N p htQ decreases in direct 
proportion to the neutral fraction of the gas, so if the IGM at z = 6 has nni/nn = 0.15 as 
suggested by Cen (2003), the quasar would have to be very highly magnified through lensing, 
in an extremely overdense region of the IGM, or extraordinarily young in order to produce 
an H II region as small as that we observe. 

What other explanations could there be for the small H II region observed in front of 
SDSS J1148+5251? One possibility is that the redshift of SDSS J1148+5251 is underesti- 
mated, which would increase the observed size of the ionized region. If the quasar redshift 
is 6.41 instead of 6.37 (Willott, McLure & Jarvis 2003), the Stromgren radius increases to 
4.7 Mpc. That would increase the N p htQ product in Eqn. (7) by a factor of 4.8 to 0.31, 
which reduces the discrepancy with the observed luminosity and expected lifetime. How- 
ever, a redshift as high as 6.41 produces a very poor match between the LBQS template 
and the observed spectrum (Fig. 3), so we consider a redshift that large unlikely. Match- 
ing the LBQS composite spectrum to the Mg II emission line detected at 2/im by Willott 
et al. (2003) reveals that z = 6.37 matches the data about as well as z = 6.41; in fact, 
from the LBQS spectrum, which does include the offsets typically seen between the redshifts 
of different emission lines, a redshift of 6.41 would appear to be an upper limit for SDSS 
J1148+5251. It will probably require higher SNR infrared spectra of this object to resolve 
the question of its true redshift. 

Another possibility is that clumpy structure in the IGM conspires to limit the size of 
the Stromgren sphere along the line of sight. Certainly there is substantial IGM structure 
near the quasar, as it can be seen in the narrow absorption lines on the blue edge of the Ly a 
emission line for wavelengths approaching the GP trough. A cloudy IGM can also modify 
the radiative transfer for ionizing radiation so that the Stromgren sphere is not fully ionized 
(Cen & Haiman 2000; Fan et al. 2002), so that gas within the ionized region contributes 
substantially to the observed Ly a optical depth. 



-27- 



However, these effects ought to be relatively small for H II regions that are still in their 
very early phase of relativistic expansion (see Appendix A). In that phase there is a surfeit 
of ionizing photons, with far more photons present than are required to ionize the neutral 
hydrogen. The expansion rate of the ionization front is limited by the light-travel time; 
recombinations are completely negligible, and even dense gas clumps are fully ionized by the 
copious radiation. 

Stopping the ionization front with a dense cloud along the line of sight is difficult for 
the same reason. For example, suppose there is a dense cloud in the IGM at z — 6.32. If 
the cloud is optically thick to Lyman continuum radiation, it will stop the Stromgren sphere 
expansion, thus producing an anomalously small ionized region. A cloud dense enough to 
stop the ionization front must have a density greater than 

(\ 1/2 / \ 1/2 , 

N ph \ _ 3 ( N ph kpc\ ( R s y 1 

A^i) = °- 12cm (lO^Tj UsMp^J ' (9) 

where «b is the case B recombination coefficient and I is the line-of-sight distance through 
the cloud. If we assume the cloud is spherical, its mass is 

M = 2xlO 6 M (J-X^fj^-Y (— . (10) 
° Vkpcy \^10 58 s J V2-8Mpcy V ; 

The cloud would have to be either dense or quite massive to stop the quasar's ionization 
front. 

Another possible explanation for the small size of the ionized region is that the edge of 
the observed emission at z — 6.32 does not indicate the edge of the ionized gas. Only a trace 
of neutral hydrogen is needed to generate a large optical depth to Ly a (eq. 2). Consequently, 
if the ionization fraction in the H II region drops below ~ 10~ 4 , the IGM becomes opaque. 
This can occur if radiation transfer in the clumpy IGM reduces the ionizing flux seen by the 
gas. This also appears unlikely in the Stromgren sphere's relativistic expansion phase, but 
more detailed calculations are required to confirm that z = 6.32 really marks the position of 
the quasar's ionization front. 

5.5. Could SDSS J1030+0524 also be lensed? 

One might wonder whether the arguments presented in favor of a foreground object 
lensing SDSS Jl 148+5251 could also be applied to SDSS J1030+0524. SDSS J1030+0524 
does show some absorption lines from intervening metal line systems (Figs. 1 and 2), but 
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none are as deep as the strongest C IV absorption line system seen in SDSS J1148+5251, 
nor do we see any evidence for associated Lya emission. The absorption spectrum of SDSS 
J1030+0524 will be discussed further by Madau & Bolte (in preparation.) 

There is no sign that continuum emission from an intervening galaxy contaminates the 
Ly« or Ly (3 GP troughs in SDSS J1030+0524. Indeed, those troughs are very black, with 
residual levels even below those detected in SDSS J 1148+5251 (see Table 2). From the flux 
limit in the Lya trough, we place a limit of wab > 27.4 on any intervening galaxy. We can 
rule out the presence of a lensing galaxy at z < 3.5 with an absolute AB magnitude brighter 
than —20 (at a rest wavelength of 2000 A). 

Finally, Haiman & Cen (2002) discussed the limits that can be placed on any lensing 
of SDSS J1030+0524 from the size of its Stromgren sphere. The quasar has an observed 
R s = 4.5 Mpc at z = 6.28, so we find (N ph /10 58 )(t Q /10 7 ) = 0.25. The ionizing photon 
luminosity is 0.13 x 10 58 s _1 or 0.61-0.87 x 10 58 s _1 (from the Elvis and Telfer templates), 
leading to age estimates of 2 x 10 7 yrs or 3-4 x 10 6 yrs. These age estimates appear plausible, 
although the higher ionizing luminosities using the Telfer et al. template do leave more room 
for lensing than was found by HC. 

5.6. Summary 

Our proposal gives a consistent picture for SDSS J1148+5251: the GP troughs of both 
Lja and Ly/5 are likely to be contaminated with light from an intervening galaxy at z ~ 5. 
The existence of this galaxy is supported by the detection of Ly a emission, C IV absorption, 
and continuum emission. There could be several such intervening systems, since both the 
emission and absorption are complex with multiple components spread over ~ 1000 km/s. 
The intervening system amplifies the quasar's light via gravitational lensing, thereby enhanc- 
ing the likelihood of discovery for the quasar. The possibility of strong lensing is supported 
by the relatively small H II region created by the quasar's ionizing radiation, which indicates 
that the quasar is probably considerably less luminous than its apparent brightness would 
indicate. 

There may be some way to reconcile the observations of SDSS J 1148+5251 presented 
here with a model based purely on IGM absorption. No such explanation is obvious to us, 
though there are so many excellent theoreticians who are interested in this problem that 
there will doubtless be many good ideas proposed shortly after the publication of these 
results! From the observational side, clearly an HST image of SDSS J1148+5251 would be 
very interesting. It could reveal the presence of multiple images (or set strong limits on their 
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absence.) An HST image taken using a narrow-band filter in the Ly/3 trough would directly 
determine whether the peaks seen at 7205 and 7236 A are extended and offset (as expected 
if they are Lya emission from an intervening galaxy) or are point-like and coincident with 
the quasar (as expected if they are the result of "leaks" in the IGM absorption.) 

If the intervening galaxy suggestion turns out to be correct, then our observations of 
SDSS J1030+0524 remain the single best measurement of the transmission of the IGM at 
z > 6. There is no light detected in the Lya or Ly (3 troughs for this object at a very low 
level. On the other hand, if SDSS J1148+5251 is shown not to have an intervening system 
contaminating its GP troughs, then the IGM toward it is in fact quite transparent, with 
both numerous high-ionization holes and a significant transmission across the whole Lya 
and Ly/3 troughs. The resolution of these questions will have to wait for the discovery of 
additional z > 6 quasars. 

Many thanks to Michael Bolte and Piero Madau for sharing their ESI exposures of 
SDSS J1030+0524 with us. We appreciate a helpful conversation with James Rhoads and 
Sangeeta Malhotra. RLW thanks Holland Ford and the Johns Hopkins University Astron- 
omy & Physics Department for their hospitality during a sabbatical, when much of this 
work was carried out, and appreciates continuing support from the Space Telescope Science 
Institute. RHB acknowledges support from NSF grant AST-00-98355 and the Institute of 
Geophysics and Planetary Physics (operated under the auspices of the U. S. Department of 
Energy by Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48). 
XF acknowledges support from the University of Arizona and an Alfred P. Sloan Research 
Fellowship. MAS acknowledges support from NSF grant AST-00-71091. 

A. Appendix: Early Evolution of a Stromgren Sphere 

The early evolution of the ionized region around a quasar is marked by a period of very 
rapid expansion, with the ionizing front moving out at nearly the speed of light. In that 
case the finite light travel time across the Stromgren sphere cannot be ignored. However, 
the size of the sphere as inferred from observations of absorption along the line-of-sight to 
the quasar turns out to have exactly the same evolution with time as one derives under the 
assumption of an infinite speed of light. This appendix derives the evolution of the size both 
as seen in the frame of the ionizing source and as seen through the line-of-sight absorption. 

We assume that the mean free path for ionizing photons in the neutral gas is very short 
and that recombinations and the expansion of the universe can be ignored. These are good 
approximations: the mean free path for 1 Rydberg photons is only 1 kpc at z = 6 and declines 
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as (1 + z) 3 . Cosmological simulations indicate that the gas clumping factor C = (n 2 H ) / (nn) 2 
is likely to be well below the value C ~ 100 required to make recombinations a significant 
factor (Gnedin k Ostriker 1997; Madau k Rees 2000; Cen k Haiman 2000). 

With these assumptions, ionizing photons propagate freely to the edge of the ionized 
region where they are immediately absorbed by the neutral gas. Figure 10 shows the relevant 
geometry. The photons that would have traveled beyond the Stromgren radius Rs (the 
shaded region) have been absorbed by neutral hydrogen atoms within Rs- The size of the 
H II region is determined by a balance between the number of photons in the outer shell and 
the number of atoms within the ionized sphere: 

^(t-^)=^L , (A1) 

where N p h is the ionizing photon luminosity and t is the age of the source. This cubic 
equation can be analytically solved to give 

^-TTJTT* • (A2) 

where 

and the time t c , which roughly marks the end of the relativistic expansion period, is 

p^y /2 =12xloV (i)" 1 (\±^y 12 (MiMiy 172 , (A4) 

\4im H c 3 y ll0 57 s-M V 7 / Wh 0.0196/ ' V ; 



where zq is the quasar redshift. Here we have explicitly included the dependence on the 
neutral fraction riui/nu and the cosmological parameters. 

The expansion law in Eqn. (A2) is shown in Figure 11. After an initial period where the 
ionization front moves out at nearly the speed of light, the expansion slows and approaches 
the standard t 1/3 law (Madau k Rees 2000, Cen k Haiman 2000). 

This is not, however, the expansion law that is observed via absorption along the line- 
of-sight. When we observe Lya absorption from the neutral gas just beyond the edge of 
the H II region, the Lya photons being detected were emitted only a short time after the 
first ionizing photons. Specifically, if the Lya photons cross the Stromgren boundary at a 
radius Rs, they were emitted when the age of the quasar was only t' — t — Rs/c (which 
is the time for light to cross the outer, shaded shell in Fig. 10). The observed Stromgren 
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Fig. 10. — Schematic diagram of early evolution of a Stromgren sphere. The outer circle of 
radius ct marks the extent (in the absence of absorption) of the ionizing radiation at a time 
t after the source turns on. The inner circle of radius Rs is the ionized Stromgren sphere. 
The number of photons in the shaded volume is equal to the number of hydrogen atoms in 
the ionized volume. 




Fig. 11. — Size versus time for early Stromgren sphere expansion. The time is scaled by t c , 
the duration of relativistic expansion (see text for details) and the radius by R c = ct c . The 
dashed line is R = ct. At late times the expansion is R ~ t 1 / 3 . 
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evolution, as inferred from absorption, is not Rs(t) versus t but rather Rs(t) versus t' . With 
some algebraic manipulation of the analytical expressions for R s and t', we can express Rs 
directly in terms of the observed (apparent) time t' as 



where R c = ct c . 

This is a surprising result (at least, it surprised us!) The expansion law for the observed 
radius is exactly the same as the expansion law derived if one completely ignores light-travel 
time effects (eq. 6). In the frame of the quasar the ionization front's expansion velocity is 
limited by the speed of light, but for the observer there is an initial period of superluminal 
expansion. The speed of light drops out of the solution along the line-of-sight because the 
delay required to allow light to travel from the source to the edge of the H II region is exactly 
compensated by the speedup that results from that edge being closer to the observer, allowing 
light originating there to reach the observer sooner than light from the quasar core. 

Note that the size of the observed H II region in Eqn. (A5) may be larger than the 
ionized region that exists at the end of the quasar's lifetime. For example, suppose a quasar 
shines for a period t c and then shuts down. According to Eqn. (A2), at t c the radius of the 
H II region is 0.7 R c . But the observer measures a radius R c (eq. A5). What is happening 
in this case is that the H II region continues to expand after the quasar shuts off until it 
reaches a maximum size R c at t — 2t c . At that point the inner sphere in Figure 10 is 
completely empty of ionizing photons and so the expansion stops. The observer is detecting 
Lya photons that were emitted just at the end of the quasar's lifetime, but by the time 
those photons cross the ionization front it has expanded to its maximum size. 

Our conclusion is that when all the light travel time effects are taken into account, 
Eqn. (A5) is the correct formula for computing the age of the quasar given the luminosity. 
Ages determined using that formula are, in fact, the age of the quasar at the time it emitted 
the photons being detected, and so such ages are direct measurements of the minimum quasar 
lifetime even during the superluminal expansion phase. 
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